
Reliable differentiation of Meyerozyma 
guilliermondii from Meyerozyma caribbica by 
internal transcribed spacer restriction 
fingerprinting 



Romi et a I. 



BioMed Central 



Romi et al. BMC Microbiology 2014, 14:52 
http://www.biomedcentral.com/1471-2180/14/52 



Romi et al. BMC Microbiology 2014, 14:52 
http://www.biomedcentral.eom/1 471 -21 80/1 4/52 



Microbiology 



METHODOLOGY ARTICLE Open Access 



Reliable differentiation of Meyerozyma 
guilliermondii from Meyerozyma caribbica by 
internal transcribed spacer restriction 
fingerprinting 

Wahengbam Romi 1,2 , Santosh Keisam 1,2 , Giasuddin Ahmed 2 and Kumaraswamy Jeyaram 1 * 
Abstract 

Background: Meyerozyma guilliermondii (anamorph Candida guilliermondii) and Meyerozyma caribbica (anamorph 
Candida fermentati) are closely related species of the genetically heterogenous M. guilliermondii complex. 
Conventional phenotypic methods frequently misidentify the species within this complex and also with other 
species of the Saccharomycotina CTG clade. Even the long-established sequencing of large subunit (LSU) rRNA gene 
remains ambiguous. We also faced similar problem during identification of yeast isolates of M. guilliermondii 
complex from indigenous bamboo shoot fermentation in North East India. There is a need for development of 
reliable and accurate identification methods for these closely related species because of their increasing importance 
as emerging infectious yeasts and associated biotechnological attributes. 

Results: We targeted the highly variable internal transcribed spacer (ITS) region (ITS1 -5.8S-ITS2) and identified seven 
restriction enzymes through in silico analysis for differentiating M. guilliermondii from M. caribbica. Fifty five isolates 
of M. guilliermondii complex which could not be delineated into species-specific taxonomic ranks by API 20 C AUX 
and LSU rRNA gene D1/D2 sequencing were subjected to ITS-restriction fragment length polymorphism (ITS-RFLP) 
analysis. Taq\ ITS-RFLP distinctly differentiated the isolates into M. guilliermondii (47 isolates) and M. caribbica (08 
isolates) with reproducible species-specific patterns similar to the in silico prediction. The reliability of this method 
was validated by ITS1-5.8S-ITS2 sequencing, mitochondrial DNA RFLP and electrophoretic karyotyping. 

Conclusions: We herein described a reliable ITS-RFLP method for distinct differentiation of frequently misidentified 
M. guilliermondii from M. caribbica. Even though in silico analysis differentiated other closely related species of 
M. guilliermondii complex from the above two species, it is yet to be confirmed by in vitro analysis using reference 
strains. This method can be used as a reliable tool for rapid and accurate identification of closely related species of 
M. guilliermondii complex and for differentiating emerging infectious yeasts of the Saccharomycotina CTG clade. 
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Background 

Meyerozyma guilliermondii is a genetically heterogenous 
complex belonging to the Saccharomycotina CTG clade 
[1]. This complex consists of phenotypically indis- 
tinguishable and closely related species namely Meyero- 
zyma guilliermondii (anamorph Candida guilliermondii), 
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Meyerozyma caribbica (anamorph Candida fermen- 
tati), Candida carpophila, Candida smithsonii, Can- 
dida athensensis, Candida elateridarum and Candida 
glucosophila [2-6]. Apart from its presence in healthy 
human [7,8], M. guilliermondii also exists in clinical 
[3,9] and environmental samples [10]. This organism 
is widely studied in various aspects due to its clinical 
importance, biotechnological applications and bio- 
logical control potential [11]. 
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C guilliermondii is regarded as an emerging infectious 
yeast of the non-albicans Candida (NAC) species group 
which accounts for 1 - 5% of nosocomial blood stream in- 
fections worldwide [9,12,13]. However, in certain geo- 
graphical regions such as Brazil, India and Italy, over 10% 
of all the candidaemia cases are caused by this species 
[14]. The threat posed by this organism is ever increasing 
due to the decreased susceptibility and emergence of 
strains resistant to antifungal drugs like polyene (ampho- 
tericin B) and azoles (fluconazole and itraconazole), lead- 
ing to mortality in candidaemia patients [9,12,15]. C. 
fermentati has been rarely found to be associated with 
candidaemia [16, 17]. But due to the poor discernability of 
C. fermentati from C. guilliermondii, they are commonly 
misidentified in clinical laboratories. Apart from being or- 
ganisms of clinical importance, M. guilliermondii and M. 
caribbica are often linked with fermented foods [18-20]. 
M. guilliermondii is known for the production of flavour 
compounds in fermented food products [21]. Further, in a 
study with soybean paste fermentation, M. guilliermondii 
and M. caribbica have been claimed for the efficient pro- 
duction of isoflavone aglycone which is a widely known 
bioactive compound for its various health promoting func- 
tions [22]. M. guilliermondii is a flavinogenic yeast which 
is known for the overproduction of vitamin B 2 (riboflavin) 
[23]. Moreover, isolates of M. guilliermondii andM carib- 
bica have exhibited great potential in the biological control 
of fungi responsible for postharvest spoilage of fruits and 
vegetables [24-26]. These yeast species with enhanced 
biological control efficacy have emerged as a potential al- 
ternative to the conventional fungicide treatment. Consid- 
ering the various importance and applications of the two 
species, there is a need for the development of accurate 
and reliable method to identify and distinctly discriminate 
the closely related species. 

Current methods of yeast identification, mostly in clin- 
ical practice, are mainly based on the conventional and 
rapidly evolving commercial phenotypic and biochemical 
methods. However, such methods are often unreliable 
for accurate identification of closely related yeast species 
[13,27]. According to recent studies, M. guilliermondii 
and M. caribbica are extremely difficult to differentiate by 
the phenotypic methods [28-31]. We also faced similar 
problem during differentiation of yeast isolates from 
soibum, an indigenous fermented bamboo shoot product 
of North East India (Additional file 1: Table SI). The 
widely used API 20 C AUX yeast identification system and 
sequencing of large subunit (LSU) rRNA gene D1/D2 do- 
main failed to give proper species-level taxonomic assign- 
ment to these isolates (Additional file 1: Tables S2 and S3). 
Moreover, the phylogenetic tree reconstructed from the 
publicly available D1/D2 sequences of different strains of 
M. guilliermondii and M. caribbica failed to discriminate 
the two species (Additional file 2: Figure SI). 



Several attempts have been made using molecular ap- 
proaches such as DNA base composition, electrophoretic 
karyotyping [6,32], multi locus sequence typing (MLST) 
[3], multi locus enzyme electrophoresis (MLEE), randomly 
amplified polymorphic DNA (RAPD) [4], sequencing of 
internal transcribed spacer (ITS) [28,30], intergenic spacer 
restriction fragment length polymorphism (IGS-RFLP) 
[29] and RFLP of housekeeping genes such as riboflavin 
synthetase gene RIBO [17] in order to resolve the misiden- 
tification. Some recent studies have claimed that the 
matrix-assisted laser desorption ionization-time of flight 
mass spectrometry (MALDI-TOF-MS) is advantageous 
over previous approaches for reliable identification of 
clinically important NAC and non-Candida yeast species 
[28,31,33,34]. Unfortunately, MALDI-TOF-MS requires 
reference spectra of accurately identified closely related 
strains otherwise the results may be erroneous. On the 
other hand, the sequence-based studies have considered 
the ITS1-5.8S-ITS2 region as universal DNA barcode for 
yeast identification [35] and the RFLP of ITS1-5.8S-ITS2 
region has successfully separated the closely related spe- 
cies in the genera Candida and Pichia [36,37] . Therefore, 
in this study, we targeted the ITS1-5.8S-ITS2 region to de- 
velop a simple RFLP method for accurate taxonomic as- 
signment of M. guilliermondii and M. caribbica. 

With this background, the aim of the present study was 
(i) to perform in silico prediction of restriction enzymes to 
discriminate M. guilliermondii and M. caribbica using the 
publicly available ITS1-5.8S-ITS2 sequences, (ii) to eva- 
luate the selected enzymes by in vitro ITS -RFLP analysis 
of ambiguously identified 55 yeast isolates for species- 
specific taxonomic assignment, and (iii) to validate the 
taxonomic assignment by ITS1-5.8S-ITS2 sequencing, 
mitochondrial DNA (mtDNA)-RFLP and pulsed field gel 
electrophoresis (PFGE) karyotyping. 

Methods 

Yeast isolates and strains 

The yeast isolates used in the present study are listed in 
Additional file 1: Table SI. These isolates were obtained 
from samples collected at different stages of indigenous 
bamboo shoot fermentation for the production of soibum 
in Manipur state of North East India [38]. The sample 
(10 g) was homogenized in 90 mL of sterile physio- 
logical saline (1 g/L bacteriological peptone, 8.5 g/L 
NaCl, pH 6.1) using Stomacher® 400 Circulator (Seward, 
Worthing, West Sussex) at 250 rpm for 3 min. The 
yeasts were isolated by serial dilution spread-plating of 
the above homogenate on yeast extract peptone dex- 
trose (YEPD) agar medium (pH 6.5) (HiMedia, Mumbai, 
India) containing 100 [ig/mL each of filter-sterilized 
ampicillin and tetracycline (Sigma- Aldrich, Bangalore, 
India), followed by incubation at 30°C for 48 - 72 h 
under aerobic conditions. All the isolates were purified 
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by sub-culturing twice on the same agar medium and 
preserved at -80°C in YEPD broth containing 10% (v/v) 
sterile glycerol (Sigma- Aldrich). For short term storage, 
the cultures were maintained at 4°C on YEPD agar. 
The type strain C. guilliermondii ATCC 6260 used for 
comparison was obtained from American Type Culture 
Collection. 

Phenotypic characterization and morphological 
observation 

Phenotypic identification of the yeast isolates was carried 
out using the API 20 C AUX yeast identification system 
(bioMerieux, New Delhi, India) following manufacturer s 
instructions. Colony and cell morphology of the isolates 
were studied using SZ-PT stereo binocular microscope 
(Olympus, Japan) and BX61 phase contrast microscope 
(Olympus). 

In silico analysis and restriction enzyme selection 

The full length ITS1-5.8S-ITS2 sequences of M. guillier- 
mondii and M. caribbica were retrieved from NCBI 
(http://www.ncbi.nlm.nih.gov/) and Centraalbureau voor 
Schimmelcultures (CBS-KNAW) yeast nucleotide data- 
bases (http://www.cbs.knaw.nl/Collections/Biolomics.aspx? 
Table = CBS + strain + database). Type strain sequences of 
the two species, C. guilliermondii ATCC 6260 [GenBank: 
AY939792.1] and AL caribbica CBS 9966 (http://www.cbs. 
knaw.nl/Collections/BioloMICS.aspx?Link = T&TargetKey = 
14682616000000137&Rec = 36291&Revert = F) were sub- 
jected to in silico PCR amplification using primers ITS1 
(5 -TCCGTAGGTGAACCTGCGG-3 ') and ITS4 (5 -TC 
CTCCGCTTATTGATATGC-3 ') [39] to trim off the un- 
targeted regions on both 5 ' and 3 ' ends of the sequences 
using the online Sequence Manipulation Suite (http://www. 
bioinformatics.org/sms2/pcr_products). Using NEBcutter, 
version 2.0 (http://tools.neb.com/NEBcutter2/), the ITS 
amplicon sequences were subjected to in silico restriction 
digestion using the commercially available type-II restric- 
tion endonucleases listed in the REBASE database (http:// 
rebase.neb.com/) [40] to select the enzymes which cut the 
two sequences differently at not more than 5 cleavage 
sites. Multiple sequence alignment of 10 additional ITS1- 
5.8S-ITS2 sequences of different strains from different 
ecological niches for each species was performed using 
Clustal X, version 2.0 (http://www.clustal.org/clustal2/) 
and BioEdit, version 7.2.0 (http://www.mbio.ncsu.edu/ 
bioedit/bioedithtml) to confirm the taxa-specificity of the 
selected restriction enzymes. 

DNA extraction 

DNA was extracted from pure cultures as cell-free DNA 
lysate using lyticase-heat lysis method. Briefly, a single 
colony of 24 - 48 h old culture from YEPD agar was in- 
oculated to 5 mL of YEPD broth supplemented with 



antibiotics, and incubated for 18 h at 30°C with shaking 
at 200 rpm. Cells were harvested from 1 mL of the cul- 
ture broth at 5,000 g for 5 min (FA-45-24-11, Centrifuge 
5424, Eppendorf, Hamburg, Germany). The cell pellet 
was washed twice with 1 mL sterile 0.5 M NaCl followed 
by sterile deionized water (Milli Q, Millipore, Molsheim, 
France). The cells were finally resuspended in 500 \iL of 
lx TE buffer (10 mM Tris-Cl, 1 mM EDTA, pH 8.0) 
containing 10 \iL of lyticase (5U/uL) (Sigma- Aldrich) 
and incubated at 37°C for 1 h. After the incubation, the 
spheroplasts were lysed by heating at 95°C for 20 min. 
The crude cell-free lysate was collected by centrifugation 
at 10,000 g for 10 min at 4°C and the DNA was quanti- 
fied spectrophotometrically (Nanodrop ND-1000, Nano- 
Drop Technologies, Inc., Rockland, USA). The cell-free 
lysate with absorbance ratio (A 2 6o/28o) of 1.8 - 2.2 was 
used for PCR analysis and stored at -20°C until 
required. 

ITS-RFLP 

ITS1-5.8S-ITS2 was amplified from the cell-free DNA 
lysate using primers ITS1 and ITS4 mentioned else- 
where. The amplification was carried out in a 25 \iL final 
reaction volume containing 50 ng of the genomic DNA 
as previously described [41]. The amplified ITS fragment 
was analyzed by 2.0% (w/v) agarose gel electrophoresis 
at 80 V in 0.5 x TBE (45 mM Tris-borate, 1 mM EDTA, 
pH 8.0) buffer to check its intactness and absence of 
non-specific amplification. The PCR product (4 |iL) was 
digested with 5 U of Taql (Promega, Madison, USA) in a 
10 (iL reaction volume at 65°C as per manufacturer s in- 
structions. The restriction patterns were analyzed by 
electrophoresis of the 10 \xL reaction volume on 2.0% 
(w/v) agarose gel in parallel with PCR 100 bp Low DNA 
ladder (Sigma- Aldrich) as molecular size standard. The 
electrophoresis was run at 80 V for 2 h in 0.5 x TBE buf- 
fer. The gel was then stained in 0.5 (ig/mL ethidium 
bromide solution for 30 min with rocking at 15 rpm 
on a platform rocker (Tarsons, Kolkata, India). After 
destaining for 30 min in sterile deionized water, the gel 
was photographed using ChemiDoc MP gel documenta- 
tion system (Bio Rad, Hercules, USA). The restriction 
fingerprints were analysed for the absence or presence 
of discriminating fragments using GelCompar II soft- 
ware, version 6.5 (Applied Maths, Sint-Martens-Latem, 
Belgium). 

mtDNA-RFLP 

A single colony of 24 - 48 h old culture from YEPD agar 
was inoculated to 5 mL of YEPD broth supplemented 
with antibiotics, and incubated for 18 h at 30°C with 
shaking at 200 rpm. The grown culture was inoculated 
into 50 mL of fresh YEPD broth (initial OD 600 = 0.1) and 
incubated in the above conditions till mid-logarithmic 
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growth phase (final OD 600 = 0.4 - 0.8). Cells of 20 OD 600 
were harvested at 1,800 g for 5 min at 4°C (A-4-81, Cen- 
trifuge 5810R, Eppendorf). The mtDNA was extracted 
as previously described [42] with some modifications. 
The cells were resuspended and washed with 5 mL of 
yeast resuspension buffer (50 mM Tris-Cl, 20 mM 
EDTA, pH 8.0) and stored at -20°C for 10 min. Lyticase 
(50 U) ( Sigma- Aldrich) was used to produce spheroplast 
and 15 \iL of 1 mg/mL RNase A solution (Sigma- 
Aldrich) was added during cell lysis. The total DNA was 
precipitated at -20°C for 1 h. After quantifying the DNA 
content spectrophotometrically, the DNA was freeze 
dried, re-dissolved in sterile deionized water to a final 
concentration of 1 |ig/|iL and stored at -20°C till further 
use. Restriction digestion was carried out on 10 \ig of 
the DNA in a 20 \xL reaction volume using 10 U each of 
Haelll and Hinfl (Promega) according to manufacturer s 
instructions. The restriction patterns were generated by 
1.0% (w/v) agarose gel electrophoresis of the 20 \iL reac- 
tion volume at 80 V in 0.5 x TBE buffer for 4 h in paral- 
lel with 1 kb DNA ladder (Promega). After staining and 
documentation, the restriction fingerprints were sub- 
jected to cluster analysis using unweighted pair group 
method with arithmetic mean (UPGMA) algorithm on 
Jaccard similarity coefficients using GelCompar II. Com- 
posite data set of the restriction digestion profiles was 
generated with 1.0% position tolerance to generate the 
clustering. Bootstrap analysis with 1,000 replicates was 
performed to indicate the branch quality. 

Electrophoretic karyotyping 

Intact chromosomal DNA for electrophoretic karyotyping 
using PFGE was prepared as previously described [32]. 
The electrophoresis was carried out in 1.0% (w/v) PFGE- 
grade agarose gel ( Sigma- Aldrich) and 0.5 x TBE buffer at 
13 - 14°C and 150 V in contour- clamped homogeneous 
electric field electrophoresis apparatus (Gene Navigator, 
Amersham Biosciences, Uppsala, Sweden). The gel was 
run for 22 h with a switch interval of 90 s for 8 h followed 
by 105 s for 6 h and finally 120 s for 8 h in parallel with 
PFGE marker (225 - 22,000 kb) from Saccharomyces cere- 
visiae strain YPH80 (Sigma- Aldrich). Staining and docu- 
mentation were performed as mentioned elsewhere. 

ITS and D1/D2 sequencing and sequence analysis 

The representative isolates from each ITS-RFLP genotype 
group were randomly selected for sequencing ITS1-5.8S- 
ITS2 and LSU rRNA gene D1/D2 domain. ITS1-5.8S-ITS2 
was PCR amplified as mentioned elsewhere. Amplification 
of D1/D2 region was carried out using primers NL1 
(5 -GCATATCAATAAGCGGAGGAAAAG-3 ') and NL4 
(5 -GGTCCGTGTTTCAAGACGG-3 ') as previously de- 
scribed [43]. The amplified products were purified using 
NucleoSpin® Extract II gel extraction kit (Machery-Nagel, 



Duren, Germany) following manufacturers instructions. 
The PCR products were sequenced using ABI 3100 Gen- 
etic Analyser (Merck, Bangalore, India) with the same 
primers used for the amplification. The sequence reads 
were validated by analysing the electropherogram data 
using ChromasLITE software, version 2.01 (http://techne- 
lysium.com.au/). To identify the closest known relatives, 
the sequences were queried with NCBI and CBS yeast nu- 
cleotide databases. The sequences obtained from both se- 
quencing and nucleotide databases were aligned using 
Clustal X algorithm and a neighbour joining tree was con- 
structed by Kimuras evolutionary distance matrix obtained 
from the multiple sequence alignment using MEGA4 
phylogenetic software. Bootstrap values for 1000 replicates 
were shown at the node of cluster branch. The sequences 
were deposited to NCBI GenBank under the following 
accession numbers: JF439366 - JF439369 and KF268351 - 
KF268354. 

Results 

in silico selection of differentiating restriction enzymes 

The full length ITS1-5.8S-ITS2 sequences of M. guillier- 
mondii and M. caribbica were retrieved from NCBI and 
CBS yeast nucleotide databases and subjected to mul- 
tiple sequence alignment followed by in silico restriction 
digestion. Three variable regions differentiating M. guil- 
liermondii from M. caribbica were identified. Seven re- 
striction enzymes (Arsl, Bfal, Bsrl, HpylSSl, HpyCW^lll, 
Mmel and Taql) which cut the variable regions differ- 
ently were identified (Figure 1A and Additional file 2: 
Figure S2). Considering the length and the number of 
polymorphic fragments with sizes greater than 100 bp 
(for easy analysis in normal agarose gel), Bfal, Mmel and 
Taql were found appropriate. Notably, commonly avail- 
able Taql gave distinct species-specific differentiation be- 
tween the two species (Additional file 1: Table S4 and 
Additional file 2: Figure S3). We also tested the selected 
three restriction enzymes (Bfal, Mmel and Taql) for dif- 
ferentiating M. guilliermondii and M. caribbica from 
other closely related members of M. guilliermondii com- 
plex (Additional file 1: Table S4). Except C. carpophila 
and M. caribbica, all other members of M. guilliermon- 
dii complex were distinctly differentiated during the 
analysis. 

Evaluation of in silico selected restriction enzymes by 
in vitro ITS-RFLP 

To validate the above in silico selection, the 55 yeast iso- 
lates of M. guilliermondii complex (which were not dif- 
ferentiated by phenotypic characterization and D1/D2 
sequencing) were analysed by ITS-RFLP using the se- 
lected Taql restriction enzyme in comparison with the 
type strain C. guilliermondii ATCC 6260. All the tested 
isolates and the type strain gave a single PCR amplicon 
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Clustal Consensus 

Saccharomyces_cerevisiae S288c (NC_001144) 

M. guilliermondii ATCC 6260 T (AY939792 . 1) 

M. guilliermondii WM10.14 (JN183444.1) 

M. guilliermondii CanR-56 (JF817286.1) 

M. guilliermondii PX-PAT (GQ497898.1) 

M. guilliermondii SD 337 (GQ334393 . 1) 

M. guilliermondii EQ (EU177579.1) 

M. guilliermondii G7A1 (DQ680842.1) 

M. guilliermondii JHSd (DQ663478.1) 

M. guilliermondii L2-8 (DQ663476.1) 

M. guilliermondii UR 9406-03 (AY939795.1) 

M. guilliermondii HJM (EF191048.1) 

M. caribbica CBS 9966 T (CBS Collection) 

M. caribbica IWBT-Y836 (JQ993381.1) 

M. caribbica CNRMA 200700593 (EU569001.1) 

M. caribbica IMUFRJ 51970 (FN428931.1) 

M. caribbica CBS 2022 (EU568913.1) 

M. caribbica WM10.15 (JN183445 . 1) 

M. caribbica UCLM 44A (GQ340430 . 1) 

Candida fermentati ATCC 22995 (AF022718.1) 

M. caribbica S58-1 (GU943488.1) 

M. caribbica H5MandK13 (FM199963.1) 

M. caribbica 1-5 (AB260139.1) 



580 590 600 610 620 630 

. . | | | | | | | | | | | . . 

★ ****★★***★★**★ *** * *** *** * * * ** 

AGTGAGTGATACTCTTTGGAGTTAACTTGAAATTGCTGGCCTTTTCATTGGATGT 

ATTGAGTGATACTCTT AGTCGGACTA GGCGTTTGCTTGAAAAGT 

ATTGAGTGATACTCTT AGTCGGACTA GGCGTTTGCTTGAAAAGT 

ATTGAGTGATACTCTT AGTCGGACTA GGCGTTTGCTTGAAAAGT 

ATTGAGTGATACTCTT AGTCGGACTA GGCGTTTGCTTGAAAAGT 

ATTGAGTGATACTCTT AGTCGGACTA GGCGTTTGCTTGAAAAGT 

ATTGAGTGATACTCTT AGTCGGACTA GGCGTTTGCTTGAAAAGT 

ATTGAGTGATACTCTT AGTCGGACTA GGCGTTTGCTTGAAAAGT 

ATTGAGTGATACTCTT AGTCGGACTA GGCGTTTGCTTGAAAAGT 

ATTGAGTGATACTCTT AGTCGGACTA GGCGTTTGCTTGAAAAGT 

ATTGAGTGATACTCTT AGTCGGACTA GGCGTTTGCTTGAAAAGT 

ATTGAGTGATACTCTT AG TCGG ACTA GGCGTTTGCTTGAAAAGT 

ATTGAGTGATACTCTT AGTCGAACTA GGCGTTTGCTTGAAAAGT 

ATTGAGTGATACTCTT AGTCGAACTA GGCGTTTGCTTGAAAAGT 

ATTGAGTGATACTCTT AGTCGAACTA GGCGTTTGCTTGAAAAGT 

ATTGAGTGATACTCTT AGTCGAACTA GGCGTTTGCTTGAAAAGT 

ATTGAGTGATACTCTT AGTCGAACTA GGCGTTTGCTTGAAAAGT 

ATTGAGTGATACTCTT AGTCGAACTA GGCGTTTGCTTGAAAAGT 

ATTGAGTGATACTCTT AGTCGAACTA GGCGTTTGCTTGAAAAGT 

ATTGAGTGATACTCTT AGTCGAACTA GGCGTTTGCTTGAAAAGT 

ATTGAGTGATACTCTT AGTCGAACTA GGCGTTTGCTTGAAAAGT 

ATTGAGTGATACTCTT AGTCGAACTA GGCGTTTGCTTGAAAAGT 

ATTGAGTGATACTCTT AGTCGAACTA GGCGTTTGCTTGAAAAGT 



B 




Figure 1 Differentiation of M. guilliermondii and M. caribbica by Taq\ digestion of ITS1-5.8S-ITS2. A: Multiple sequence alignment of 
representative ITS1 -5.8S-ITS2 sequences of the two species obtained from NCBI GenBank and CBS yeast database showing position of Taql 
recognition site (highlighted) which distinctly differentiated the two species. B: Taql restriction digestion profile of ITS1-5.8S-ITS2 amplicons 
obtained from some of the representative isolates. Lane 1: C. guilliermondii ATCC 6260; Lane 2-12: isolates of M. guilliermondii genotype group 
MG (A1S10Y1, A2S10Y1, A3S9Y1, A2S9Y1, A3S11Y1, A3S2Y1, A3S6Y1, A2S6Y1, A1S9Y1, Kw3S3Y1 and Kw2S11Y2); Lane 13-20: isolates of 
M. caribbica genotype group MC (A1S10Y2a, A1S10Y3, A1S10Y5, Kw3S2Y1, Kw2S3Y1, Kw3S3Y3, Kw3S3Y4 and Kw1S7Y2); Lane M: PCR 100 bp 
Low DNA ladder (Sigma-Aldrich). 



of molecular size of 607 bp. As predicted by the in silico 
analysis, Taql ITS-RFLP distinctly differentiated the iso- 
lates into two genotype groups. Forty seven isolates pro- 
duced M. guilliermondii-specific pattern (MG), while the 
remaining eight isolates generated M. caribbica-specific 
pattern (MC) (Table 1). Examples of Taql ITS-RFLP 
profiles differentiating the above two species are shown 
in Figure IB. 

Validation of ITS-RFLP assay 

The above ITS-RFLP based discrimination of M. guillier- 
mondii and M. caribbica was further confirmed by 
ITS1-5.8S-ITS2 sequencing, mtDNA-RFLP fingerprint- 
ing and PFGE karyotyping (Table 1). The ITS sequences 
of the isolates in each genotype group MG and MC 



matched with the sequences of the type strains C. guil- 
liermondii ATCC 6260 and M. caribbica CBS 9966 with 
99.6% and 99.8% similarity respectively. The sequences 
between the two groups were 99% identical showing 
only 5 nucleotide differences which were the same as 
shown by the above type strain sequences. Unlike Dl/ 
D2 region, the ITS sequences formed distinct cluster of 
M. guilliermondii and M. caribbica during phylogenetic 
analysis (Figure 2). The ITS sequences of M. guillier- 
mondii strains PX-PAT, CanR-56 and SD 337; M. carib- 
bica strains UCLM 44A and IWBT-Y836 showed 1 - 3 
nucleotide differences which was reflected as divergence 
within the clusters. The mean evolutionary divergence of 
0.0131 between the two clusters was 6 times more than 
the divergence within each cluster. 
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Table 1 Differentiation of ambiguous 55 yeast isolates obtained from soibum into Meyerozyma guilliermondii and 
Meyerozyma caribbica 



Group Representative Taxonomic designation 



(Number 
of isolates) 


strains 


API 20 C 

AUX* 


Taq\-\TS- 
RFLP 


Sequencing 
LSU D1/D2 ITS1-5.8S-ITS2 


mtDNA- 
RFLP 


Karyotyping 


MG (47) 


A1S10Y1, 


M. 


M. 


M. guilliermondii/ / M. caribbica 


M. guilliermondii 


M. 


M. 




Kw2S11Y2 


guilliermondii 


guilliermondii 


(JF439368, JF439369) + 


(KF268351, KF268352) 


guilliermondii 


guilliermondii 


MC (08) 


Kw1S7Y2, 


M. 


M. caribbica 


M. guilliermondii/ / M. caribbica 


M. caribbica (KF268353, 


M. caribbica 


M. caribbica 




Kw3S2Y1 


guilliermondii 




(JF439366, JF439367) 


KF268354) 






Type strain 


ATCC 6260 


M. 


M. 


M. guilliermondii/ / M. caribbica 


M. guilliermondii 


M. 


M. 






guilliermondii 


guilliermondii 


(AJ508562.1) 


(AY939792.1) 


guilliermondii 


guilliermondii 



*No identification data for M. caribbica is included in the database. 
+ GenBank accession numbers. 



The mtDNA-RFLP using Haelll and Hinfi distinctly seg- 
regated the yeast isolates into M. guilliermondii and M. 
caribbica. mtDNA-RFLP profile-based dendrogram formed 
two clusters (Figure 3) similar to the ITS-RFLP groups. 
Between the two enzymes used, Hinfi showed higher poly- 
morphism than Haelll. Electrophoretic karyotyping also 
distinctly discriminated the above two species (Figure 4). 
The species-specific mtDNA-RFLP pattern suggested that 
the isolates of each group belonged to only one strain 
(Figure 3). Whereas electrophoretic karyotyping brought 
out strain level diversity in both the groups which con- 
firmed that multiple strains of M. guilliermondii and 



M. caribbica were involved in the indigenous bamboo shoot 
fermentation (Figure 4 and Additional file 2: Figure S4). 

Discussion 

In recent times, the frequency of emerging infectious 
diseases caused by the opportunistic yeast species of 
NAC and non-Candida groups has increased in im- 
munosuppressed patients [12,44]. This is linked with the 
indiscriminate use of broad-spectrum antifungal drugs 
and global climate change [45-47]. Most of these closely 
related yeast species are often misidentified by the 



M. guilliermondii ATCC 6260 1 (AY939792. 1 ) 

■M. guilliermondii PX-PAT(GQ497898.1) 
M. guilliermondii WM 1 0. 1 4 (JN 1 83444. 1 ) 
M. guilliermondii HJM (EF 1 9 1 048. 1 ) 
M. guilliermondii UR 9406-03 (AY939795.1) 
M. guilliermondii JHSd (DQ663478.1) 
-M. guilliermondii CanR-56 (JF8 1 7286. 1 ) 
■M. guilliermondii SD 337 (GQ334393. 1 ) 
Kw2SllY2 (KF268352) 
M. guilliermondii EQ (EU 1 77579. 1 ) 
A1S10Y1 (KF268351) 
M. guilliermondii G7A1 (DQ680842.1) 
M guilliermondii L2-8 (DQ663476. 1 ) 
M. caribbica S58-1 (GU943488.1) 
M. caribbica CNRMA 200700593 (EU569001 . 1 ) 
M. caribbica H5MandK 1 3 (FM 1 99963. 1 ) 
KwlS7Y2(KF268353) 
M. caribbica IMUFRJ 5 1 970 (FN42893 1.1) 
M. caribbica 1-5 (AB260 1 39. 1 ) 
M. caribbica CBS 2022 (EU568913.1 ) 
Kw3S2Yl (KF268354) 
C.fermentati ATCC 22995 (AF022718.1) 
J M. caribbica WM 1 0. 1 5 (JN 1 83445. 1 ) 
M caribbica CBS9966 1 ' 
W. caribbica UCLM44A (GQ340430.1) 
M caribbica I WBT-Y836 (JQ99338 1.1) 
S.cerevisiae S288c (NC_001144) 



M. guilliermondii 
(MG) 



M. caribbica 
(MC) 



Figure 2 Neighbour-joining (NJ) phylogenetic tree showing taxa-specific separation of M. guilliermondii from M. caribbica. The tree was 

constructed based on the evolutionary distance calculated using Kimura-2 parameter from the nucleotide sequence of ITS1-5.8S-ITS2. The 

percentage of replicate trees in which the associated taxa clustered together in the bootstrap test (1000 replicates) is shown next to the branches 

for values >40%. The bar represents 1% sequence divergence. GenBank accession numbers are mentioned within the parentheses. S. cerevisiae 

was the outgroup in the analysis. T =Type strain. 
\ J 
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i, 





Figure 3 mtDNA-RFLP based dendrogram showing distinct clustering of M. guilliermondii and M. caribbica. The dendrogram was 
constructed using UPGMA algorithm on Jaccard similarity coefficients generated from Hae\\\ and Hinf\ restriction digestion profile of mtDNA of 
some of the representative isolates. Value at each branch node indicates the branch quality with 1000 bootstrap replications. The scale represents 
the similarity. 



conventional phenotypic, biochemical and antibiotic sus- 
ceptibility methods. Thus, accurate identification of 
these species has become essential for clinical manage- 
ment and effective antifungal therapy. 

In the present study, we showed that the developed ITS- 
RFLP method was reliable and consistent for distinct dif- 
ferentiation of closely related AL guilliermondii from AL 
caribbica for which phenotypic methods and D1/D2 se- 
quencing were inconclusive. Our results also indicated 
that sequencing of both D1/D2 and ITS regions will in- 
crease the resolution of species identification which can be 
further improved by multigene sequence-based phylo- 
genetic approach [3,48,49]. However, the presence of 
incorrectly identified, insufficiently annotated and non- 
updated entries in the public nucleotide databases may 



M 



10 11 



Figure 4 PFGE karyotype patterns of isolates belonging to M. 

guilliermondii and M. caribbica genotype groups. Lane 1 : C. 
guilliermondii ATCC 6260; Lane 2 - 3: M. guilliermondii isolates 
A1S10Y1 and Kw2S1 1Y2; Lane 4 - 1 1: M. caribbica isolates A1S10Y2a, 
A1S10Y3, A1S10Y5, Kw3S2Y1, Kw2S3Y1, Kw3S3Y3, Kw3S3Y4 and 
Kw1S7Y2; Lane M: 5. cerevisiae PFGE marker (Sigma-Aldrich). Right 
arrow indicates the co-migrating chromosomal doublets showing 
strain level diversity. 



underestimate the resolving power of these taxonomic 
markers [50] . Out of the 29 sequences of LSU rRNA gene 
for M. guilliermondii available in NCBI GenBank, 17 
sequences (58.62%) clustered with M. caribbica type 
strain CBS 9966 [GenBank: EU348786] (Additional file 2: 
Figure SI). 

The choice of appropriate restriction endonucleases is 
critical for RFLP experiments. The commonly used Cfol, 
Haelll and Hinfl [37,41] failed to segregate M. guillier- 
mondii from other species of the same genus during in 
silico and in vitro ITS -RFLP analysis. Our results indi- 
cated that in silico selection of restriction enzymes using 
the publicly available sequences from various strains of 
the target species is a better approach than randomly 
selecting the previously described and commonly used 
enzymes. This approach has been proven to be highly ef- 
fective and reproducible [36,51-53], and many online re- 
sources have been developed for this purpose [54-57]. 
Clinical isolates of Candida famata and Candida pal- 
mioleophila were also frequently misidentified as AL 
guilliermondii [30,31]. In silico analysis confirmed that 
the developed ITS -RFLP method can also discriminate 
these species (data not shown). This in silico selection 
approach can be effectively applied to other cryptic yeast 
species of clinical importance for the development of 
RFLP based diagnostic tools. 

The developed method of ITS -RFLP using Taql differ- 
entiated M. guilliermondii and M. caribbica at species 
level. This method is simple, rapid and reliable in com- 
parison to the commonly used sequencing methods. The 
entire analysis starting from DNA extraction to ITS- 
RFLP profiling can be completed within 8 h. Further 
studies using higher number of strains of these two 
species from different clinical sources are required to 
confirm the robustness of this method for diagnostic 
applications. Though the combination of ITS -RFLP 
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D1/D2 domain. Table S4. List of the selected type-ll restriction 
endonucleases that differentiated M. guilliermondii from M. caribbica and 
other species of M. guilliermondii complex during in silico restriction 
digestion of the ITS1 -5.8S-ITS2 amplicon sequences. 

Additional file 2: Figure SI. Neighbour-joining phylogenetic tree 
based on LSU rRNA gene D1/D2 sequences showing taxa-nonspecific 
segregation of M. guilliermondii strains. The tree was constructed based 
on the evolutionary distance calculated using Kimura-2 parameter from 
the representative nucleotide sequences of M. guilliermondii and M. 
caribbica (position 13 to 308 of LSU rRNA gene of 5. cerevisiae CBS 1 171, 
GenBank Accession No. AY048154.1). The percentage of replicate trees in 
which the associated taxa clustered together in the bootstrap test (1000 
replicates) is shown next to the branches. The bar represents 1% 
sequence divergence. GenBank accession numbers are mentioned within 
the parentheses. 5. cerevisiae was the outgroup in the analysis. T =Type 
strain. Figure S2. In silico identified restriction enzymes which distinctly 
differentiated M. guilliermondii from M. caribbica. Multiple sequence 
alignment of representative ITS1 -5.8S-ITS2 sequences of various strains of 
the two species obtained from NCBI GenBank and CBS yeast database 
showing position of identified Ars\ (A), Bfa\ (B), fori (C), Hpy188l (D), 
HpyCH4\\\ (E), and Mme\ (F) restriction recognition sites (highlighted) 
which distinctly differentiated the two species. The nucleotide position 
was based on the sequence of the in silico PCR amplicon of ITS1-5.8S- 
ITS2 of 5. cerevisiae strain S288c (NC_001 144) including gaps generated 
during multiple sequence alignment. C. fermentati is the anamorph of 
M. caribbica. T =Type strain. Figure S3. In silico restriction digestion 
profile of M guilliermondii and M. caribbica ITS1 -5.8S-ITS2 amplicon. The 
theoretical restriction digestion profile was generated using NEBcutter, 
version 2.0 (http://tools.neb.com/NEBcutter2/). Lane G: M. guilliermondii 
ATCC 6260; Lane C: M. caribbica CBS 9966; Lane M: 100 bp DNA ladder. 
Figure S4. Strain level diversity of M. guilliermondii revealed by PFGE 
karyotyping. Lane 1-13: Isolates A3S2Y1, Kw1S2Y1, Kw3S3Y1, A3S6Y1, 
A2S6Y1, A1S9Y1, A1S9Y5, A2S9Y1, A3S9Y1, A3S9Y9, A2S10Y1, A2S10Y4 
and A3S1 1Y1. White arrow indicates the polymorphic chromosomal 
band. 



profiles generated by TaqI, Bfal and Mmel differentiated 
other closely related species of the M. guilliermondii 
complex from M. guilliermondii and M. caribbica during 
in silico analysis, it is yet to be confirmed through 
in vitro analysis using reference strains. 

In our present study the use of either Haelll or Hinfi for 
mtDNA-RFLP was sufficient to differentiate the two spe- 
cies, but failed to give strain level differentiation even 
when both enzymes were used. This demonstrated that 
mtDNA-RFLP can also be used for distinct differentiation 
of closely related species. The Hinfi mtDNA-RFLP pattern 
of our M. guilliermondii isolates was similar with the 
mtDNA restriction pattern of M. guilliermondii strains 
isolated from wineries in Alentejo, Portugal [58]. This 
genotype was linked with the production of flavour com- 
pound, 4-ethylphenol in wine. The major phenolic flavour 
compound (4-methylphenol) detected from fermented 
bamboo shoot product, soibum (Singh NR: unpublished 
observations) might also have originated from M. guillier- 
mondii. Future study is required to characterize the 
flavour compound producing strain for starter culture de- 
velopment. Though fresh bamboo shoots are highly per- 
ishable, the fermented bamboo shoot can be preserved up 
to one year after fermentation without any deterioration 
or change in its organoleptic character. This long term 
preservation may be linked with the dominant presence of 
M. guilliermondii which has been reported as an efficient 
biological control agent [24,25]. Being an emerging infec- 
tious yeast, the presence of M. guilliermondii in fermented 
food is a great concern regarding the safety of its con- 
sumption. Further study in strain level is required to 
unravel the pathogenic potential of M. guilliermondii asso- 
ciated with soibum fermentation. 

Conclusions 

In this study, we described an ITS-RFLP method devel- 
oped through an integrated approach of in silico selection 
of restriction enzymes and in vitro validations for distinct 
differentiation of frequently misidentified M. guilliermon- 
dii from M. caribbica, which can be used as an alternative 
or an adjunct to ITS sequencing. This method may be 
used for rapid and accurate identification of emerging in- 
fectious yeasts of the Saccharomycotina CTG clade. This 
approach can also be used for other closely related species 
complex when phenotypic methods and D1/D2 sequen- 
cing are ambiguous. 

Additional files 



Additional file 1: Table SI. List of the 55 yeast isolates used in the 
present study. Table S2. Carbon substrate assimilation pattern of 
representative strains of M. guilliermondii complex using API 20 C AUX 
yeast identification system. Table S3. Taxonomic assignment of isolates 
belonging to M. guilliermondii complex by sequencing of LSU rRNA gene 
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